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Abstract 


The structure of chiral symmetry restorations at finite temperature is 
thoroughly investigated in the supersymmetric Nambu-Jona-Lasinio model 
with a soft supersymmetry breaking term. It is found that the broken chiral 
symmetry at vanishing temperature is restored at sufficiently high tempera¬ 
ture in two patterns, i. e., the first order and second order phase transition 
depending on the choice of the coupling constant G and the supersymmetry 
breaking parameter A. The critical curves expressing the phase boundaries 
in the G — A plane are completely determined and the dynamically generated 

fermion mass is calculated as a function of temperature. 

PACS numbers: ll.10.Wx, 11.15.Pg, 11.30.Qc, 11.30.Rd, 12.60.Jv 


Typeset using REVTpX 


1 


The supersymmetry is an essential ingredient of unified field theories for elementary 
particles so that it is worth trying to study a supersymmetric version of composite Higgs 
models such as Technicolor model. We shall consider a supersymmetric composite model 
in the present communication with an interest in the early universe and take into account 
finite temperature effects. As a prototype of the composite model we pick up the Nambu- 
Jona-Lasinio (NJL) model which is known to be useful to investigate the mechanism of 
the dynamical symmetry breaking [|l]. The supersymmetric version of the NJL model is 
easily constructed. Unfortunately, however, in the supersymmetric NJL model the chiral 
symmetry is strongly protected to keep the boson-fermion supersymmetry (SUSY) and hence 
the dynamical chiral symmetry breaking does not take place 0. If a soft SUSY breaking term 
is added to the SUSY NJL Lagrangian, the dynamical breakdown of the chiral symmetry 
is brought about for sufficiently large SUSY breaking parameter A [Q. The reason for this 
is simple: The large A implies the large effective mass of the scalar components of the 
superfields so that quantum effects due to the scalar components get suppressed compared 
with that of the spinor components. Thus the model becomes closer to the original NJL 
model which allows the dynamical fermion mass generation. Stating the same substance in 
a different way we realize that the boson by acquiring its mass term forces supersymmetry 
to make balance so that the fermion mass is generated dynamically. 

The SUSY NJL model with soft SUSY breaking is useful to study the mechanism of 
the dynamical chiral symmetry breaking within the framework of supersymmetry. If we 
take the model seriously as a prototype of the unified field theory, it is natural to extend 
the argument to take into account circumstances of the finite temperature and spacetime- 
curvature as in the early universe. It should, however, be noted that SUSY is not a good 
symmetry for finite temperature |4|]. This fact may discourage studies of SUSY models 
at finite temperature. In our argument we consider the supersymmetry of the model in 
the strict sense only at vanishing temperature and focus our attention on the dynamical 
breaking of the chiral symmetry in the same model at finite temperature. The breaking of 
SUSY at finite temperature is kinematical since it is caused essencially by the difference of 
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the statistical distribution functions for bosons and fermions. 


A natural SUSY extension of the NJL model is characterized by the following Lagrangian 

0 . 

f - - G 

Csnjl = J d 4 0 [$+$ + + 

-AW 2 (<5 + $ + + $_$_)], (1) 

where d A 6 = d 2 6d 2 6. The chiral superfield is represented by component scalar (4>±) 
and spinor (ijj±) fields and auxiliary field F± so that $±( 1 /) = <j>±{y) + V^0^±(y) + F±(y)9 2 
with y = x + iOa^Q. The chiral superfields $ + and <f>_ belong to the multiplet N and N in 
SU(N) respectively. Note that a soft SUSY breaking term is introduced in Eq. ([[]) following 
Buchmuellcr and Ellwanger [Q. 

In the following arguments it is convenient to introduce auxiliary chiral superfields H 
and S and use a Lagrangian equivalent to Eq. (i 

C = J d 4 0[($+$ + + - A 2 e 2 e 2 ) + — HH\ 

r N 

+ [J d 2 e (-HS - S'$+$_) + /i.c.]. (2) 

Written in terms of component fields the Lagrangian takes the following form where we have 
kept only relevant terms for calculating our effective potential in the leading order of the 
1/A expansion: 

C = —</>+( □ + \4>s\ 2 + 2k 2 )4>+ — </>-(□ + |<As| 2 + 2k 2 )0_ 

-cr + i7 5 7r)T - ^|0 S | 2 . (3) 

The field efts is related to the ordinary scalar (a) and pseudoscalar (7r) held such that <f>s = 
a + Z7r, and T = ^-) T . 

We would like to calculate the effective potential for Lagrangian (||) in the leading order of 
the 1/N expansion at finite temperature. We adopt the Matsubara formalism |5|] to introduce 
temperature effects in our calculation of the effective potential. We simply replace the 
integral in energy variable by the summation where we apply the periodic boundary condition 
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for boson fields and the anti-periodic boundary condition for fermion fields respectively. 
The calculation of the effective potential is essentially the same as the one described in our 
previous paper | 6 j. The summation can be easily performed and after some algebra the 
following result is obtained, 

^(|0s| 2 ) = ho(|0s| 2 ) + Vg(\(fis\ 2 ), (4) 

where Vo(| 0 s| 2 ) is the effective potential for vanishing temperature, 

VMM 2 ) = ^ + 2/AU±( 

— \Jk 2 + \(j> s \ 2 — Vk 2 + A 2 + Vk 2 ), (5) 


which agrees with the one obtained by performing the ko integration after Wick rotation 
in the effective potential given at vanishing temperature, and Vg(| 0 s| 2 ) is the temperature 
dependent part of the effective potential that is given by 


VMM 2 ) = ~j 


d D ~ l k .l + e _/ V fc 2 +l*sl 2 

ln(- 


PJ ( 2 tt) d - 1 “ v 1 4 - e -/Vfc 2 

1 _ p -/Vfc 2 +A 2 

), 


X- 


( 6 ) 


l _ e -/?V fc2 +|0sl 2 +A2 

where D is the space-time dimension (Although D = 4, we keep D to be arbitrary in Eqs. 
(|j) and @ for later convenience), (3 is the inverse of the temperature T in unit of kp, — \ 
with ks the Boltzmann constant. Note that the effective potential (|J) is normalized so that 


E(0) = 0. 

ft is easy to see that the effective potential (|D reduces to Vo(| 0 s| 2 ) in the zero tempera¬ 
ture limit f3 —*■ 00 . If we let A —> 0 in the zero temperature effective potential (||), we are left 
only with the tree term as is expected by the boson-fermion cancellation in the SUSY model. 
In the SUSY limit A —> 0 the finite temperature effective potential (pj) still retains the term 
of quantum corrections which is composed of the bosonic and fermionic statistical distribu¬ 
tion functions. This result reflects the well-known fact that supersymmetry is not a good 
symmetry at finite temperature |3j. At sufficiently low temperature both the Fermi-Dirac 
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and Bose-Einstein distribution function 1 /(exp (3yjk 2 + \(ps\ 2 ± 1) are well approximated by 
the Maxwell-Boltzmann distribution function exp {—(3yjk 2 + |</>s| 2 ) and hence the quantum 
correction term in Eq. (jl|) disappears. Thus supersymmetry is a good symmetry not only 
at zero temperature but also in the low temperature region. 

Although the phase structure of our model at finite temperature is studied by direct 
numerical estimates of the effective potential (j4|), it is more convenient to deal with the gap 
equation V'{\(j)s\ 2 ) — 0 where the prime denotes the differentiation with respect to \4>s\ 2 '■ 

Hl<fe| 2 ) = dV/d\<f> s \ 2 , 

1 f d D ~ l k ( 1 1 . 

G J ( 2 vr )^— l( ^2 + | 0 s |2 ^ 2 + | 0 s |2 + A 2 ) 

^ f d D ~ X k r 1 1 

+ J (2ti) D - 1 y/fc 2 + \(j) S \ 2 e /V fc2 + lti?| 2 + 1 

-7 1 - , 1 -] = 0. (7) 

y]k 2 + |0 5 | 2 + A 2 e V fc2 +l7sl 2 +A 2 _ x 

In order to see the situation as clearly as possible we start with the discussion in the case 
D = 3 where the gap equation (^) after performing the integration is much simpler than the 
one for D = 4 and there is no divergence in the integration in Eq. (|?|). The gap equation 
for D = 3 is obtained in an analytical form by performing the integration in Eq. (j?]) : 

27r /- /- 2 1 + e _/3 V , l7sl 2 

<8) 

By numerical inspections we find that Eq. (|8|) allows at most two nontrivial solutions for 
1 4>s | 2 depending on the choice of parameters G, A and (3. Thus the effective potential (£|) 
with D = 3 may have zero, one and two extremum (extrema) respectively. Hence, when the 
temperature 1/(3 is increased, we will have possibilities of experiencing no phase transition, 
the 2nd order phase transition and the 1st order phase transition. 

The effective potential V(\((>s\ 2 ) for D = 3 is calculated analytically such that 

l^lfel 2 ) = ^ + A 3 - (Ifel 2 + A 2 ) 3 / 2 ] 

+ ^s(I ( / ) s| 2 )d=3) (9) 
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where we have suppressed writing explicitly the temperature dependent term (the last term 
on the right hand side of Eq. ([])) which is rather lengthy and can be expressed by using 
the polylogarithms. 

At vanishing temperature the chiral symmetry in the model is broken J| if sufficiently 
large is the parameter A appearing in the soft SUSY breaking term. In fact by observing 
the behavior of the gap equation V^(\(j)s\ 2 ) = 0 we can confirm the statement. In three 
dimensions for vanishing temperature [3 —> oo the last term on the left hand side of Eq. (||) 
disappears and we see that only one nontrivial solution for \(j)s\ 2 exists when 2n/G < A. 
The curve on the G — A plane given by the equation 


GA/(2vr) = 1 


( 10 ) 


divides the whole G — A plane into two regions which represent the unbroken and broken 
chiral symmetry respectively. 

If the parameters G and A are kept in the unbroken region, we remain in the unbroken 
region when temperature is increased. If the parameters G and A are in the broken region, we 
experience the chiral symmetry restoration as temperature gets high enough. This situation 
is easily confirmed by observing the gap equation (^]) and the effective potential (0. In the 
present situation there are two types of the chiral symmetry restoration : the first-order and 
second-order phase transition. These two cases can be distinguished from each other by the 
conditions U'(0) = 0, U"(0) < 0 and U'(0) = 0, U"(0) > 0 respectively. Hence the condition 
that = 0 be a point of inflection for the effective potential, U'(0) = U"(0) = 0, 

divides these two cases and gives a curve on the G — A plane. The curve is the critical 
curve on the G — A plane which divides the region of the broken phase into the regions of 


the first and second order phase transition. Written in parameters G and A the condition 
U'(0) = U"(0) = 0 reads 


GA f3 A 1 — e~ /3A ’ 

(3A/2 = coth(/3A/2). 


2t r 2,2 

— = 1 ——- In- 


( 11 ) 


( 12 ) 
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Here (3 A/2 is determined by Eq. (Jl2|) and is equal to 1.19968. The critical curve distin¬ 
guishing the second order phase transition from the first order one is given by the equation 
obtained by substituting the above value for /3A/2 in Eq. (|TT|), i.e., 

GA/(2vr) = 2.91631. (13) 

The whole G — A plane is divided by two hyperbolas (jlOj) and (fl~3|) into three regions which 
correspond to the unbroken, broken (2nd order) and broken (1st order) region respectively 
as shown in Fig. 1. Note here that as can be seen in Eq. (j9|) the effective potential for D = 3 
reduces to that of the ordinary NJL model by letting A —► oo. The physical reason for this 
result is simple: In the limit of large A the scalar components of the superfield carry a large 
effective mass and decouple from the spinor components leaving only the contribution of 
the NJL fermions. It should, however, be remarked that the gap equation (|8|) diverges as 
A —> oo. This result corresponds exactly to the fact that the gap equation for the ordinary 
NJL model diverges linearly for D = 3. In fact it is given by 

§ = A-v/Sf-|Ml+e-' , v / ^F) (14) 

and corresponds to Eq. (|3|) with large A replaced by Lambda. In the ordinary NJL model 
for three dimensions the bare coupling constant G is to be kept small for large A. This fact 
precisely reflects the property that G gets smaller as A becomes larger along the critical 
curve. 

We then go over to the case of the four dimensional space-time. For D = 4 the situation 
is almost the same as in D = 3. The main difference comes from the presence of the 
divergence for D = 4 in the integration of Eqs.(£|) and (0). We regularize the divergence 
by introducing the three dimensional cut-off A in the momentum integration. The effective 
potential for D = 4 at vanishing temperature Vo (|0s| 2 ) is calculated analytically while the 
temperature dependent part of the effective potential Vg(|0,s| 2 ) is estimated numerically. In 
these calculations all the parameters in the model are scaled by the cut-off A. By direct 
numerical estimates we find that the effective potential (j4|) with D = 4 represents the chiral 
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symmetry restoration as temperature gets high enough. In Fig. 2 (a) and (b) we see that 
the chiral symmetry restoration is of the 2nd as well as 1st order phase transition depending 
on the choice of the parameters G and A (In our previous paper the existance of the 2nd 
order phase transition was overlooked 0). 

The gap equation <91//<9|0s| 2 = 0 reads 


1 rz - rz - z , 1 + \/l + X + (5 

- yjl + x + vl + x + 8 — (x + 8) In - - 

2 a y/x + 8 


+x In 


1 + y/1 + x An 2 dVp(\(t)s\ 2 )D=A 


= 0, 


(15) 


ix A 2 d\4>s\ 2 

where a, 8 and x are defined by a = GA 2 /(87r 2 ), 6 = A 2 /A 2 ,x = |</>s| 2 /A 2 . For (3 —► oo the 
last term on the left hand side in Eq. (0) drops out leaving only temperature independent 
terms and we have the gap equation at zero temperature. It is not difficult to show that 
Eq. ( |T5|) without the temperature dependent term allows a nontrivial solution for x when 
the parameters a and 8 satisfy the condition, 


< l-VTTS + Sln^^^. (16) 

2a v 8 

Thus the dynamical fermion mass is generated and the chiral symmetry is broken dynam¬ 
ically at vanishing temperature if the above inequality is satisfied. The boundary of the 
region given by Eq. ( 0 ) in the a — 8 plane is the critical curve dividing the symmetric and 
broken phase as shown in Fig. 3. 

Just in parallel with the case of D = 3 we find three possibilities in which the gap equation 
(0) allows no solution, one nontrivial solution and two nontrivial solutions respectively 
depending on the choice of parameters G, A and (3. It is easily found by observing the 
behaviors of the gap equation and the effective potential that there are two types of phase 
transitions, the first order and second order, as temperature increases. As in the case of 
D — 3 the critical curve distinguishiding the region of the 1st order phase transition from 
the region of the 2nd order phase transition is given by solving simultaneous equations for 
G and A corresponding to the following equations, 


y'(o) = y"(o) = 0. 


( 17 ) 













The critical curve derived from Eq. (|T7|) is shown by a dotted line in Fig. 3. 

The full and dotted line divide the whole G — A region into three regions of the unbroken 
chiral symmetry, the broken chiral symmetry with the 2nd order phase transition and the 
broken chiral symmetry with the 1st order phase transition respectively. The dynamical 
fermion mass generated by the above chiral symmetry breaking is obtained by numerically 
solving the gap equation (|T5D as is shown in Fig. 4. 

One of the direct physical consequences of our approach may be found in the held of elec- 
troweak baryogenesis in early Universe ||. In this connection it is interesting to see whether 
composite Higgs masses can be taken large enough to account for the present experimental 
observation and to construct realistic composite Higgs models for describing the baryogen¬ 
esis 0. As a physicsl application of our approach we may consider the supersymmetric E e 
GUT model [Hj in which the symmetry breaking patterns can be precisely investigated by 
using our results. It is also interesting to note whether our results presented in the paper 
is subject to any change if we introduce the contribution of order 1/N or higher. Such 
calculation has been done in Ref. [ITT] in the case of the nonsupersymmetric NJL model in 
three dimensions and an extension of their result to the case of the SUSY NJL model should 
be straightforward. The investigation in this direction will be left for future works. 

After the completion of our work we became aware of a report discussing the similar 


analysis as in the present paper |12|. The conclusion of Ref. |12], however, is quite different 
from ours. 

The authors would like to thank K. Kikkawa and C. S. Lim for useful comments on 
supersymmetry at finite temperature and T. Inagaki, S. Mukaigawa and M. Tanabashi for 
enlightening discussions. 
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FIGURES 



FIG. 1. Critical curves on G — A plane in D = 3: ’S’ represents the region of the parameters G 
and A in which the chiral symmetry is not broken at any temperature, ’II’ represents the region of 
the parameters in which the broken chiral symmetry at vanishing temperature is restored at high 
temperature as the second order phase transition and ’I’ represents the region of the parameters 
in which the broken chiral symmetry at vanishing temperature is restored at high temperature as 
the first order phase transition. 
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FIG. 2. The (3 dependence of the effective potential (a) for a = 25 and 6 


0.1 and (b) for 


a = 90 and <5 = 0.1. 
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FIG. 3. Critical curves on G — A plane in D = 4: The explanation of symbols ’S’, ’II’ and ’I’ 
is the same as the one in Fig. 1. 
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FIG. 4. The dynamical fermion mass as a function of /3 (a) for a 


25 and 5 


0.1 and (b) for 


a = 90 and <5 = 0.1. 
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